Fractional diffusion limit for collisional kinetic 
equations: A moments method 



This paper is devoted to hydrodynamic limits of linear kinetic equa- 
tions. We consider situations in which the thermodynamical equilibrium 
is described by a heavy-tail distribution function rather than a maxwellian 
distribution. A similar problem was addressed in [14] using Fourier trans- 
form and it was shown that the long time/small mean free path behavior 
of the solution of the kinetic equation is described by a fractional dif- 
fusion equation. In this paper, we propose a different method to obtain 
similar results. This method is somewhat reminiscent of the so-called "mo- 
ments method" which plays an important role in kinetic theory. This new 
method allows us to consider space dependent collision operators (which 
could not be treated in [14]). We believe that it also provides the relevant 
tool to address nonlinear problems. 
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1 Introduction 

Our goal is to study the asymptotic behavior as e goes to zero of the solution 
f^{x^v^t) of the following kinetic equation: 
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when the operator L is a Unear relaxation collision operator of the form 

Hf) - / [a{x,v,v')f{v')-a{x,v',v)fiv)]dv' (2) 

(0(e), which is such that linig^o ^'(e) = will be chosen later). As usual, the 
collision operator is decomposed into a "gain" term and a "loss" term as follows: 

L(/) = K{.f) - vf 

with 

K{f){x,v) = / a{x,v,v')f{x,v')dv' 
and the collision frequency v defined by 

v{XtV) ^ I (t{x,v' ,v) dv' . 

This singular perturbation problem is very classical. The underlying goal 
is the derivation of macroscopic models that describe the evolution of a cloud 
of particles (represented, at the microscopic level, by the distribution function 
f^{x,v,t)) for small Knudscn number (of order e) and large time (of order 
0{e)~'^). The collision operator (2) is one of the simplest operator that mod- 
els diffusive, mass-preserving interactions of the particles with the surrounding 
medium. We recall that under reasonable assumptions on the collision kernel 
a{x,v,v') (see [6]), one can show that there exists a unique cquilibriimr function 
F{x,v) > satisfying 

L{F) = and / F(.t, v) dv = 1 a.e. x £ M^. 

JR« 

Formally, Equation (1) then leads to 

lim f e kcr(i) = {p{x, t)F{x, v) ; p : M.^ x (0, oo) > M}. 

Our goal is to show that the density p{x, t) is solution to an equation of hydro- 
dynamic type. 

The derivation of hydrodynamic limits for kinetic equations such as (1) was 
first investigated by E. Wigner [18], A. Bensoussan, J.L. Lions and G. Papan- 
icolaou [2] and E.W. Larsen and J.B. Keller [13] and it has been the topic of 
many papers since (see in particular C. Bardos, R. Santos and R. Sentis [1] 
and P. Degond, T. Goudon and F. Poupaud [6] and references therein). In [6], 
P. Degond, T. Goudon and F. Poupaud considered very general collision oper- 
ators of the form (2). When F decreases "quickly enough" for large values of 
\v\ (it is often assumed that is a Maxwellian distribution function of the form 
F{v) = Cexp( ^^^'^ )), they proved in particular that for 6{e) = e^, f'^{x,v,t) 
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converges, when e goes to zero, to a function of the form p{x,t)F{x,v) where 
the density p{x,t) solves a drift-diffusion equation: 

dtp-div,{D\7,p + Up) = 0. (3) 

This resuh is proved in [6] under some assumptions on F that guarantee in 
particular that the diffusion matrix D and the coefficient [/, which depends on 
F and a, arc finite. 

When is a power tail (or heavy tail) distribution function, satisfying 

^(") ~ j^;^ a« 1^1 °« (4) 

for some a > 0, the diffusion matrix D in (3) might however be infinite. In 
that case, the diffusion limit leading to (3) breaks down, which means that the 
choice of time scale 8(e) ~ e'^ was inappropriate. It is the goal of this paper to 
investigate such situations. 

Power tail distribution functions arise in various contexts, such as astrophys- 
ical plasmas (see [17], [15]) or in the study of granular media (see [9], [5] for 
the so-called "inelastic Maxwell model" introduced in [3]). We refer to [14] for 
further references concerning the relevance of power tail distribution functions 
in various applications. 

In [14], S. Mischler, C. Mouhot and the author addressed this problem in 
the space homogeneous case (that is with a independent of x). It was shown 
that when F satisfies (4) and the collision frequency v satisfies 

v{v) i^ol'^l'^ as \v\ oo, 

then for some values of a and /? (a > and /3 < min(a, 2 — a)), the choice of 
an appropriate time scale 9{e) leads to a fractional diffusion equation instead 
of (3). Let us give a precise statement in the simplest case which corresponds 
to collision kernels satisfying 

i^iF{v) < ct{v,v') < V2F{v) 

(i.e. /3 ~ 0) when (4) holds with a G (0, 2). Then, taking 6{e) ~ e", it is shown 
in [14] that the fimction f^{x,v,t) converges to p{x,t)F{v) where p{x,t) solves 
the following fractional diffusion equation: 

dtp + K{-A)''/^p = 0. (5) 



We recall that the fractional power of the Laplacian appearing in (5) can be 
defined using the Fourier transform by 

^((-A)"/V)(fc) = ikmrm- 

Alternatively, we have the following singular integral representation: 

(-A)"/V(X)^«„PV/^^|M_^.„. (6) 
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We refer the reader to Landkof [12] and Stein [16] for a discussion of the prop- 
erties of fractional operators and singular integrals (we only need to know the 
definitions above for the purpose of this paper) . 

Anomalous diffusion limits for kinetic models was first investigated in the 
case of a gas confined between two plates, when the distance between the plates 
goes to (see [4], [10], [7], [8]). In that case, the Hmiting equation is still a 
standard diffusion equation, but the time scale is anomalous {0(e) ~ ln(£^^)) 
and the particles travelling in directions nearly parallel to the plates are respon- 
sible for the anomalous scaling. A fractional diffusion equation such as (5) was 
obtained for the first time as a diffusive limit from a linear phonon-Boltzmann 
equation simultaneously by A. Mellet, S. Mischlcr and C. Mouhot in [14] and 
by M. Jara, T. Komorowski and S. OUa in [11] (via a very different probability 
approach). 

The method developed in [14] to establish the result quoted above (and more) 
relies on the use of Fourier transform with respect to x, and an explicit com- 
putation of the symbol of the asymptotic operator. Unfortunately, this method 
is rather complicated to implement when the collision operator depends on the 
space variable (the Fourier transform of L{f^) involves some convolutions), and 
virtually impossible to use in a nonlinear framework. In this paper, we thus 
propose a different method, which is closer to the so-called "moment method" 
classically used to study hydrodynamic limits of kinetic equations. This method 
relies on the introduction of an appropriate auxiliary problem and a weak for- 
mulation of (1). More precisely, the corner stone of the proof is to multiply 
Equation (1) by a test function x'^ix,v,t), solution of the auxiliary equation 

i^{x, v)x'' - £v ■ V^x" = i^ix, v)(p{x, t) 
where 1^9 is a smooth test function. 

This method allows us to consider space dependent collision kernels (which 
could not be treated previously) and we believe that it will provide the relevant 
tools to address non-linear problems (such as equations coupled to Poisson's 
equation via an electric field or equations involving non-linear collision opera- 
tors). 

In the next section, we present the main results of this paper, starting with 
a simple study case which we will use to present the method in a simpler frame- 
work (Theorem 2.1). We then address the case of general collision operators 
(Theorem 2.2). Section 3 and 4 are devoted to the proof of the main theorems. 



2 Main results 
2.1 The simplest case 

The simplest framework corresponds to collision kernels (t(x, v, v') of the form 

a{v,v') = F{v) 
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with F : (0, oo) positive function normalized so that 

/ F{v) dv = 1. 
In that case, the colhsion operator reduces to: 

/ f{v')dv'F{v)-fiv), (7) 



and we have 

ker(L) = {p{x,t)F{v) ; p : x R ^ M}. 
We will also assume that the equilibrium distribution function F{v) satisfies 

FeL°°{R^), and = i^(-iO for all w G (8) 

(note that the symmetry assumption is always crucial in deriving diffusion equa- 
tions from kinetic models). 

In that case, classical arguments (see for instance [6]) show that if we choose 
9{e) — in (1), then /" converges as e goes to zero to p{x,t)F{v) with p 
solution of: 

dtp - dxXi^ijdx.p) = 

where 

Kij = / ViVjF{v)dv. 

However, when F satisfies (4) with a S (0, 2), we get ku = oo, so this limit does 
not make sense. Our first result addresses this case: 

Theorem 2.1. Assume that L is given by (7) with the normalized distribution 
function F{v) > satisfying (8) and 



F{v) — > Ko > as \v\ oo (9) 



for some a G (0,2). Let f'^{x,v,t) be a solution of (1) with 9(e) ~ e" and 
fo e ly^iR^ X R^) nLi(K^ X R^), fo > 0. 

Then f- converges weakly in L°°{0,T; Lj,_^{R^ x R^)) to p{x,t)F{v) with 
p solution of 

1 dtp+K{-Ar/^p)^o ^^^^ 

\ p(x,0) = po{x) 



where 



CN,a Jo 



dz 



and pq = Jjjjv /o(a;, z;) (the constant cjv.q is the constant appearing in the 
definition of the fractional Laplace operator (6)). 
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In this result (and in the rest of this paper), x M.^) denotes the 

weighted set equipped with the norm: 

Note that Theorem 2.1 is proved in [14] using Fourier transform. We give 
a complete proof of this result in Section 3, because it provides the simplest 
framework to present our new "moment" method. In the next section, we 
obtain new results using this same method. We point out that the fractional 
Laplace operator will be obtained here in its singular integral form, rather than 
its Fourier symbol. This explains that the diffusion coefficient k is given by very 
different formulas here and in [14]. 



2.2 General space-dependent collision operators 

We now consider general collision operators of the form (2). We are mainly 
interested in situations in which the cross-section a depends on x (this is the 
case that could not be treated in [14]). Assuming that the equilibrium function 
F is still independent of x, the classical diffusion limit, corresponding to the 
time scale 9{e) ~ e^, leads to 

dtp - diva:{D{x)Vxp) = 

where 

-^i^) — / V (gi xix,v)F{v) dv, with x solution of L*{x) = 
Jr" 

{L* denotes the adjoint operator to L in i^(R^)). As before, it is readily 
seen that for some equilibrium distribution function F and collision kernel a, 
the matrix D{x) may be infinite. These are the situations that we wish to 
investigate. 

Before stating the result, we need to make the conditions on a and F precise. 
The first assumptions are very standard: 

Assumptions (Al) The cross-section a{x,v,v') is non negative and locally 
integrahle on M?^ for all X . The collision frequency v{x,v) — I a{x,v' ,v) dv' 
satisfies 

v{x, -v) = iy{x, v) >0 for all x, v eR^ x . 

Assumptions (A2) There exists a function F{v) G L^(R^) independent of x 
such that 

L{F) = 0. (11) 



6 



Furthermore, the function F is symmetric, positive and normalized to 1 ; 
F{-v) = F{v) > for all v eR^ and / F{v) dv = 1. 

Note that under classical assumptions on ct, the existence of an equilibrium 
function is in fact a consequence of Krein-Rutman's theorem (see [6] for details) . 
In general, however, this function may depend on x while we assume here that F 
depends only on v (we discuss in Section 2.3 the case of a;-dependent equilibrium 
function). A particular case in which these assumptions are satisfied is when a 
is such that 

Vu,w' e a{x,v,v') = b{x,v,v')F{v), with b{x,v',v) = b{x,v,v'). (12) 

In that case, we say that a satisfies a detailed balanced principle or a micro- 
reversibility principle, while the more general assumption (11) is called a general 
balanced principle. 

The next assumptions concern the behavior of F and v for large \v\. In order 
to keep things simple, we will make strong assumptions on v and F that makes 
the limiting process easier (see Remark 2.3 below): 

Assumptions (Bl) There exists a G (0, 2) and a constant kq > such that 
|w|"+^F(v) — > Ko as \v\ oo. (13) 

Assumptions (B2) There exists vi and V2 positive constants such that 

viF{v) < a{x,v,v') < V2F{v). 

In particular, integrating this condition with respect to v, we deduce: 

Q<vi< v{x, v') < V2, for all x, u' G x M^. 

Furthermore, we assume that there exists a function vo{x) (satisfying v\ < 
i^o{^) ^1^2) such that 

i/(x,z;) — > i^a{x) as \v\ 00. (14) 

uniformly with respect to x. Finally, we assume that v is with respect to x 
and 

||-Da;!^(a;,w)||L~(R2«) < C. 

We are now ready to state our main result: 
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Theorem 2.2. Assume that (Al), (A2), (Bl) and (B2) hold and let f''{x,v,t) 

-2 

/o >0 



6e a solution of (1) with 9{e) = and /o G i|,-i(M^ x R^) n ^1(10^ - 
i 0. 

Then f converges weakly in L°°(0, T; L^_i (M^ x M^)) to p{x,t)F{v) with 
p solution of 

j dtp + KoC{p) = 
1 p{x,0) = pa{x) 



(15) 



where pa = /j^n fo{x,v)dv and C is an elliptic operator of order a defined by 
the singular integral: 

with ^ 

i-oo —z / I'Q {{1 ~ s)x + sy) ds 
l{x,y) ^VQ{x)vQ{y) j -'o dz. 

Jo 

In view of (B2), it is readily seen that there exist 71 and 72 such that 

< 71 < 7(2;, ?/) < 72 < 00. 

In particular the operator C has the same order as the fractional Laplace op- 
erator (— A)"/^. It is the fractional equivalent of the divergence form elliptic 
operator — div (DV • ) which is typical of conservation laws. Note also that C is 
self-adjoint since 7(0;, j/) — j{y, x). 

Remarks 2.3. 

(i) As in [14], we could replace Assumption (Bl) by F{v) = F(j{v)£{\v\) where 
Fq satisfies (Bl) and £{\v\) is a slowly varying function (such as a logarithm, 
power of a logarithm or iterated logarithm). This would complicate the proof 
slightly, but can be handled as in [14]. In particular, as in [I4], the slowly 
varying function £ does not affect the asymptotic equation, but the time scale 
9{e) would have to be adjusted (9(e) = £"£(e~^) instead of 9(e) = e"j. 

(ii) Assumption (B2) is far from necessary, and as in [I4] macroscopic limits can 
be derived with degenerate collision frequency satisfying v{x,v) ^ vq{x)\v\^ as 
\v\ 00. In that case the operator C depends both on the asymptotic behavior 
of F and that of a. A tedious but straightforward adaptation of the proof of 
Theorem 2.2 would yield a similar result with a' = instead of oi (see [I4J). 

In particular, we must take 9{e) ~ and the operator C would be replaced by 



1-/3 Jr« \x- yY^+°' 
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2.3 Space dependent equilibrium distribution functions 
and acceleration field 

In the results above, we have assumed that the equihbrium distribution func- 
tion F was independent of the space variable x (even though the cross-section 
may depend on x). However, one might want to consider equilibrium functions 
depending on x and satisfying 

\v\°'+^F{x,v) — > Koix) as \v\ oo (16) 

instead of (13) with a e (0,2) and ko{x) > in L°°(R^). Another natural 
question concerns the effects of an external acceleration field E{x,t). In that 
case, it turns out that the correct scaling is given by 

J e°'dtr + ev- + e°'~^E -VyP ^ L{P) x G R^, n e R^, <> 
\P{x,v,{)) = h(x,v) xeR^, t;GR^ ^^"^^ 

and formally, the method described in this paper easily shows that /^(x,w,i) 
converges to p(x, t)F{x, v) with p solution of 

r dtp + CM+AiY{Ep) = Q ^^^^ 
\ P(a;,0) = Pq{x) 

where C is the elliptic operator defined in Theorem 2.2 and po = /rjv fo{x^ v) dv. 

However, in that case, the derivation of the crucial estimates (29) and 

(30) fails. Deriving the appropriate estimates (in LlogL) to rigorously justify 
this limit will be the goal of a forthcoming paper. 



3 Proof of Theorem 2.1 

In this section, we give a detailed proof of Theorem 2.1. We recall that the 
collision operator is given by (7) where F[v) > satisfies the normalization 
condition / F[v) dv = 1 and (8). 

This proof is organized as follows: We first recall the derivation of the re- 
quired a priori estimates. Then, wc introduce the auxiliary problem (23) and 
establish the main properties of the solution of that problem. Step 3 is both the 
simplest and most interesting step: Using the solution of the auxiliary problem 
in the weak formulation of (1) we show how the asymptotic problem arises as 
£ ^ 0, at least formally. The final step is devoted to the rigorous justification 
of the limits needed in Step 3. 



Step 1: A priori estimates for First of all, the maximum principle 

and the conservation of mass (consequence of the fact that J L{f) dv = for all 
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/) yield: 

||/''(i)llLi(R«xR") = / r{x,v,t)dxdv = \\fo\\Li(jsiN yj^N -J for alH > 0. 

Next, we want to prove that converges to a function of the form p{x, t)F{v). 
The foUowing computations are very classical and we only recall them here 
for the sake of completeness: We write = p'^F + with = /jg^ dv. 
Multiplying (1) by f^F^^ and integrating with respect to x and w, we get the 
following equalities: 

e"^ / ^^^F-Uvdx = / LiDf'F-Uvdx 



WfF^{rfF-^]dvdx 
/ [f - p" Ff F~^ dvdx, 



from which we deduce the two estimates 

\2 



and 



sup/ illihj)Ldvdx< [ J±dvdx = \\fo\\l2 , (19) 
/ / [f-p^FrF-Uvdxdt<-\\M\l. . (20) 



We can thus write 

r{x,v,t)=p'{x,t)F{v)+g'{x,v,t), 

with 

||5llL^ ^ <Ce"/2. (21) 
Cauchy-Schwarz inequality also gives: 



p^t,x)= I ^F'/^-dv<{ I '-i^dv 



so that 



sup/ p-(t,.)'dx< ll/olli. . (22) 

t>0 JR« -F"^ 



We deduce that converges weakly in L°°(0, T; L2,_i(R^ x M^)) to a func- 
tion p{x,t)F{x,v) where p is the weak limit of in L°°(0, T; i^(M^)). 
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step 2: The auxiliary problem. This is the key step of the proof: For 

a test function ^ £ P(R^ x [0,oo)), we introduce x^{x,v,t) G L^((0, oo) x 
solution of 

-evV,,x' = ^{x,t). (23) 

This equation is actually easy to integrate, and wc check that x^(a;, v, t) is given 
by the explicit formula: 



X''{x,v,t)= / e ''ipix + evz,t)dz. (24) 
Jo 

The function is thus smooth and bounded in L°°. Furthermore, we have: 

poo 

\x^-(p\ = / e~^[ip{x + evz,t) - ip{x,t)] dz 
Jo 

< \\Dip\\L^e\v\, 

and thus 

x'^ixjVjt) — > (p{x,t) as e — > 

uniformly with respect to x and t. However, this convergence is not uniform 
with respect to v, so wc will need the following lemma: 

Lemma 3.1. Let ip e V{R^ x [0, oo)), and define x" by (24). Then 

/ F{v)[x^ {x,v,t) — ip{x,t)\ dv — > uniformly w.r.t. x and t. 

and 

I F{v) [9tX^(a;, w, t) ~ dt(p{x, i)] dv > uniformly w.r.t. x and t. 

Furthermore, 

IIX^IlL|,(R2iVx(0,oo)) < lkllL2((o,oo)xR") 

and 

l|c?tX'^llL|,(R2«x(0,oo)) < l|c^t'PllL2((o^oo)xR«)- 

The space i|.((0, oo) x M^^) denotes the weighted space with weight 
Fiv). 

Proof of Lemma 3.1. First of all, we write 

/ F{v)[x'^{x,v,t)-Lp{x,t)\dv = I F{v) / e~''\ip{x+£vz,t)-Lp{x,t)\dzdv. 
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For some large M, we then decompose this integral as follows: 

/•OO 

F{v) / e^^[(p{x + evz,t) — (p{x,t)] dz dv 
Jo 

<\\Dip\\L^ [ F{v) I e-''£\v\zdz 

J\v\<M Jo 



dv 



+2||^||loo / F{v) / e-^'dzdv 

J\v\>M Jo 

<\\D^\\L^eM + 2\MLo. [ F{v)dv. 

J\v\>M 

Since F is intcgrable with respect to u, it is readily seen that for any (5 > 0, we 
can choose M such that /|„|>^,/ ^(^) < ^ ^'^d then choose e so that eM < 6 and 
deduce that for e small enough we have 



< C6 



where the constant C only depends on ip. The first limit in Lemma 3.1 follows. 
The other limit can be proved similarly (note that t is only a parameter here) . 

In order to prove the L^p estimates, we note that dz = 1 and so 

/>oo 

Ix^f < / e^^ipix + evz^t)"^ dz. 
Jo 

We deduce: 



IX^Ili|,(R2"x(0,oo)) 



/ / F\x'\^dxdvdt 
Jo Jr2« 



< I F e~''ip{x + evz,tydzdxdvdt 



< 



Fiv) / e ''\ML^R'^x{o.oo))dzdv 

J 



^ II'i^IIl2(rNx(0,oo))- 

The last inequality is proved similarly. □ 

Step 3: Weak formulation and formal passage to the limit. Multi- 
plying (1) by x^ and integrating with respect to x, v, t, we get: 



dtx'^ dx dv dt — e" / fo{x,v)x'^ {x,v,0) dx dv 







p'Fx' - rx' + rev ■ V,x' dx dv dt 
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which, using (23), yields: 

f^dtx^ dx dv dt — I fo{x,v)x'^ix,v,0)dxdv 







/ p^Fx" - r'p(x,t)dxdvdt 

/ p"^ Fx'^ dx dv dt — / p^ip{x,t)dxdt. 

J 



We deduce: 



f^dtx'^dxdvdt— I fo{x,v)x'^{x,v,0)dxdv 

/ [x'^(x, w, i) — iy9(x, i)] dvdxdt 

^R" JR" 



or 



dtx^ dx dv dt — I fo{x,v)x'^ {x,v,0) dx dv 

JR2« 



/ / p''C'{(p)dxdt (25) 

Jo JR" 



C'icp) = / F(«) [x%x, V, t) - ifiix, t)] dv. 



with 



The rest of the proof of Theorem 2.1 consists in passing to the hmit e — *■ in 
(25). FormaUy, we immediately see that the left hand side should converge to 



pdt^pdxdt— j po{x)ip{x,0) dx 

JR" JR" 

(this will actually be a consequence of Lemma 3.1). Passing to the limit in 
the right hand side of (25) is the most interesting part of the proof since the 
nonlocal operator should now appear in the limit of C^: Using formula (24) for 
X^ we get: 

C'{(p)=e-°'[ [ e-''F{v)[ip{x + evz,t)~ip{x,t)]dzdv. 
The change of variable w = vez and (9) then yields (formally) 



1 



£^(^) = / / j-—e-''-F{w/{£z))[ip{x + w,t)-ip{x,t)]dzdw 



R« 



(e^) 

1 



1 _ iez)^^"- 

~^ ^ I |jV+a ['^(^ + W,t)~ V?(X, t)] dz dw 



|7V+Q 



[(p{x + w,t) ~ y^i^i t)] dz dw 



R" JO 

k(-A)"/V- (26) 
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So we should expect the right hand side of (25) to converge to 




and we recover the weak formulation of (10) as the limit of (25). 

Of course, those arc just formal computations, so the next and final step 
consists in rigorously justifying those limits. Note that we will have to be par- 
ticularly careful with the right hand side, since the limit of C^{(p) involves a 
singular integral which can only be defined as a principal value (see Proposi- 
tion 3.2). 



Step 4: Rigorous passage to the limit in (25). First, we justify the 
limit in the left hand side of (25). We note that 

f^dtx^ dx dv dt 

FOtx"^ dx dv dt - / / g'^dtx'^ dx dv dt. 



/O JR2« Jo JK2JV 

The second term is obviously bounded by 

llff''llL^_i((0,oo)xR2«)||i9tX^||L|((0,oo)xR2") < Ce"^^ 

(where we use (21) and Lemma 3.1). The first term can be written as 

p''ix,t) / F{v)dtx'' dvdxdt 

JR" JR« 

which converges, using Lemma 3.1, the weak convergence of in L°°{Q, oo; L^(R^)) 
and the fact that pe is bounded in L°°(0, cxd; L^(K")), to 

\ p{x,t) / F{v)dt'p dv dx dt ^ I I p{x,t)dt'p dx dt. 
Jr" Jr" Jo Jr" 

Proceeding similarly with the initial data term, we deduce that the left hand 
side of (25) converges to 



pdttpdxdt— / po{x)ip{x,0) dx. 
Jr" Jr™ 

In order to pass to the limit in the right hand side of (25), we need the following 
lemma, which is the rigorous justification of (26): 

Proposition 3.2. Assume that the conditions of Theorem 2.1 hold and that x'^ 
is defined by (24). Then 

C'iip) / Fiv) [x^(x, v, t) - ^{x, t)] dv 
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converges as e goes to zero to 



-k(-A) / (^) = KCTV.a^'V^ / 17 ^|jv+^ d-V 



\y-x 



wher 



Ko 



\z\"e-^dz. 



C-N.a JO 

Furthermore, the convergence is uniform with respect to x and t. 

Proposition 3.2, the weak convergence of p'^ and the fact that is bounded 
in L°°{0, oo; L^(M")) allow us to pass to the limit in the right hand side of (25). 
We deduce: 



pdtfdxdt— / po{x)ip{x,0) dx =: — / p — —{~A)"^'^ipdxdt 



which is the weak formulation of (10). This completes the proof of Theorem 
2.1, and it only remains to prove Proposition 3.2. 



Proof of Proposition 3.2. We fix C > (to be chosen later). Then, using (24), 
we write: 



C'{^) = / F{v)[x'[x,v,t)-^{x,t)]dv 

e~^i^(u) Yp{x + ewz, t) — (p{x, t)] dz dv 



R" Jo 

£ " / / e~^ F{v)[ip{x + evz,t) ~ (p{x,t)] dz dv 
\v\<c Jq 



v\>C Jo 



+e-- I I e 

\v\>C Jo 

h+h + h 



Ko 



\N+a 



F{v) 



[(p{x + evz, t) — (p{x, t)] dz dv 



Ko 



\N+a 



[ip{x + evz, t) — fix, t)] dz dv 



Using the fact that F{—v) = F{v), we can write 
l^il = , , 

J|ti|<C Jo 

/ e~^ Ifix + svz, t) — ip{x, t) — ezv ■ Wip{x, t)\ dz dv 
\v\<c Jo 



e " / / s F{v)[if{x + £vz,t) — Lp{x,t) - ezv ■VLp{x,t)\dz dv 
J\v\<c Jo 



< lli^llL-£" 



< \\F\\L^\\D'if\\L^e-'' / / e-'\ezv\'dzdv 

J\v\<C Jo 



< \\F\\L^\\D'if\\L^C'e' 



(27) 
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which goes to zero as e goes to zero since a < 2. 
Next, a simple change of variable w = ezv gives 

h = r I e-^\ez\''+--^[ip{x + w,t)~ip{x)]j^dwdz 

Jo J\w\>cez |w| ^ \ez\ 

= '^0 / 1^1 e / LTpv+S '^'^dz. (28) 

JQ J\w\>Cez 1^1 

We note that by the very definition of Cauchy Principal Value, we have, for 
every z > 0, 

^.^ ^{x + w,t) -^{x,t) ^ ^{x + w,t) -^{x,t) 



= -(-A)"/V- 

However, this limit is not obviously uniform with respect to z or .t, so we have 
to work a little bit more in order to pass to the limit in (28). 

In order to prove the convergence of I2 , we recall that the fractional Laplace 
operator can also classically be written as 

-(-A^/V = PV/ [^{x + w)-v{x)]y^ 

[ip{x + w) — (p{x)]- 



[ip{x + w) - ifiix) - W(p{x) ■ w] 

,w\<l \W\ 

where all the integrals are now defined in the usual sense (no principal values). 
Using a similar decomposition in (28), we get: 

f^/'^^'^ _ f dw 
I2 = Ko \A°'e i Vp{x + w,t) - Lp{x,t)\ .j^,^ dz 

Jo J\w\>l I'^l 

/■^/(^"' _ /■ dw 
+Ko \z\°-e I [ip{x + w,t) ~ ip{x, t) - Vip{x) ■ w] dz 

Jo JCez<\w\<l 1^1 

+Ko kre-M [ip{x + w,t)-ip{x,t)]-—^dz 

Jl/iCe) J\w\>Cez \w\^+'^ 

(note that we need to split the z integral to account for the case when Cez > 1). 
The last term can be bounded by 

m\L^ \zre-^dz 
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which goes to zero as e goes to zero (uniformly with respect to x and t), while 
the first two terms converge, uniformly with respect to x and respectively to 

Ko kre""^ / [(p{x + w,t)~(p{x,t)] .j^.^ dz 

Jo "'|ti'l>l 1^1 

= KCN,a[ Mx + W,t)- lf{x, t)] I J^j^^^ 



w\>l \W\ 



and 



Ko I kre"^ / [ip{x + w,t)- ip{x, t) - Vip{x) ■ w] ^ dz. 
J|ui|<i 1"^! 

= n CN,a / + w,t) — ip{x, t) — V(^(a;) • lu 



(Note in particular that the integrand in the second term is bounded by '^^^+a-r ) 
which is integrable at w ~ Q). 
Wc deduce: 

lim/2(a;,i) = -n{-l^)°'/'^(p{x,t) 

e— >0 

and the convergence is uniform with respect to x and t. 

Finally, wc show that I3, can be made as small as we want: For any i5 > 0, 
we can choose C large enough so that 



< , for aU \v\ > C. 



(note that C was arbitrary up to now) . Proceeding in the same way as with I2 , 
we can thus show: 

limsup|/3| <5|C(^)|- 
Since this holds for any (5 > 0, the proposition follows. □ 



4 Proof of Theorem 2.2 

The proof of Theorem 2.2 is very similar to that of Theorem 2.1, so we will 
only present in detail the steps that are significantly different. We recall that 
f solves (1) with 9{e) = e". 

We start with recalling the following classical a priori estimates (see Ap- 
pendix A for the proof): 

Lemma 4.1. The solution of (1) is hounded in L°° (0, oo, L^(M^^)) and 
L°°(0, oo; L^_i (M^^)) uniformly with respect to e. Furthermore, it satisfies: 

r = p'F{v)+g', 
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where the density p*^ = J^j^ dv and the function are such that 

IIp^IIl~(o,cx),l2(r«)) < ||/o||l|,_i (29) 

and 

II/IIl^(o,oo;L^_,(k-) < C \\fo\\Ll_, (30) 

In particular converges L°°(0,T; L^)-weak to p, and converges L°°{0,T; L^_i)- 
weak to f ~ p{x,t)F{v). 

Next, we introduce the auxiliary problem corresponding to the general col- 
lision operator: For a test function if G P(R^ x [0, oo)), we define x^(a;,w,i:) 

by 

fix, w)x^ - ev ■ W^x' = i^ix, v)(p{x, t). (31) 

This equation is slightly more complicated than (23) because of the x-depcndence 
of the collision frequency i'{x^v). However, introducing x(s,x,w,t) = xt^^ + 
sev, x,v,t), a simple computation yields the following formula: 



^{x + evs, v) ds 

x'^{x,v,t)= j e -Jo ly/^x + evz,v)(p{x + evz,t) dz (32) 



(it is relatively easy to check that this function indeed solves (31)). We note 
that 

/•oo — / ^{x + evs,v) ds roo 

e Jo v{x + evz,v)dz = I e~'^ du = 1 (33) 

Jo Jo 

and so x"(x,w,t) — > ip as £ goes to zero. Furthermore, we can prove the 
following equivalent of Lemma 3.1 (we recall that is bounded from above and 
from below by positive constants): 

Lemma 4.2. For any ip e P(R^ x [0, oo)), i/x^ is defined by (32), then 
F{v)\^X^{x,v,t) — </j(.T,t)] dv — > uniformly w.r.t. x and t. 
and 

F{v)[dtx'^ {x,v,t) ~ dt(p{x,t)^ dv > uniformly w.r.t. x and t. 

Furthermore, 

Ix''IIl|,((0,Oo)xR2«) < C'||(^||i2((o,C5o)xR") 

and 

l|f^iX^llL|,((0,oo)xR2«) < C||9t<p||i2((o^oo)xRW)- 
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We now have to derive the equivalent of the weak formulation (25). For that 
purpose, we multiply (1) by and integrate with respect to x,v,t. We get: 



f^dtx^dxdvdt — e°' I fQX^{x,v,Q)dxdv 
- / / KiDx' -lyfx' + rev-V^x'dxdvdt 

Jo JR2" 

which, using (31) and the fact that K{F) = vF, yields: 



dtx'^ dx dv dt ~ I fox'^{x,v,0)dxdv 
fK*{x') - fv(p{x,t)dxdvdt 

JE2iv 

f [K*{x') - K*{ip)] dxdvdt 



A'(/^) [x^ - (fi] dxdvdt 
Jm.^" 



Kip'F) [x' - ifi] dx dv dt + / K{g^) [x' - ip] dx dv dt 

JK^n J 

p^vF [x^ ~ If] dx dv dt + / K{g^) [x" - dx dv dt (34) 
' J 

The proof of Theorem 2.2 now consists in passing to the limit in (34). As in 
the proof of Theorem 2.1, the left hand side is relatively easy to handle, with 
the help of Lemma 4.2. It converges to 



pdt(pdxdt~ / poip{x,0)dx. 

Passing to the limit in the right hand side of (34) is thus once again the most 
interesting part of the proof. We note that there is an additional term involving 
g"^ which did not appear in (25). This term will be shown to converge to zero 
as stated in the following lemma: 

Lemma 4.3. For any test function (p G T>{M.^ x [0, oo)), let x*^ be defined by 
(32). Then 

hm I I K{g') [x" - 'p] dx dv dt = 

Finally, the key proposition, which is the equivalent of Proposition 3.2 is the 
following: 

Proposition 4.4. For any test function Lp G T>{M.^ x [0, oo)), let x^ be defined 
by (32). Then 

hme-" [.ix,v)Fivnx^ix,v,t)-^M] dv^-Ci,) (35) 
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where C is defined in Theorem 2.2. Furthermore, this limit is uniform with 
respect to x and t. 

We leave it to the reader to check that Lemma 4.3 and Proposition 4.4 allow 
us to pass to the limit in (34) and thus complete the proof of Theorem 2.2. We 
now turn to the proofs of these two results: 

Proof of Lemma 4.3. First of all, we write (using Assumption (B2)): 
\K{g)\ < f a{v,v'Mv')\dv' 



< i^2F{v) ( / F{v) dv' 



1/2 



\-^dX' 

Fiv') 



< i^2J^||g||L^_i(MiV). 

for all X and v. Hence 

/ / K{g'')[x'' - if] dxdvdt 

Jo JK21V 



< 1^2 



^R" 



115 IIlI 



Fix" -(p\ dv] dxdt 



^ '^2||5'^||l^_j((0,oo)xR2«) 



"'R™ 



/ F|x^-</?|dw) dxdt] 



1/2 



(36) 



with 



and 



Since 



"^1 = / / (/ F\x'~ip\dv] dxdt 

Jo ^K" \J\v\<e-^ J 

^2 = / / I / F\X' dv] dxdt. 

Jo JK" \"'|t'|>£^l / 



Jo 



X + evs, v) ds 

v{x + evz, v)\ip{x + evz, t) — ip{x, t)\ dz 
oo — / v{x + evs, v) ds 



e "'o 



^{x + £uz, f )|£uz| dz 



< C\\Dip\\L^\evl 
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we get (using the fact that F{v) < i \n+c )■ 



Ji < CWDifW^ / Fiv)\ev\dv 

W|t,l<e-i 



2 



To estimate J2, we first note that (33) impUes: 



Jo Jk" 



00 r roQ 



i'{x + evs, v) ds 

< / e Jo ^{x + evz,v)\ip{x + evz,t)\'^ dz dx dt 



JW^ Jo 

00 



poo 

< / e~^l^zy2||¥'lli2((o,oo)xR")'^^ 

Jo 



10 

- ^II'^IIl2((o,oo)xI1 



and so, using the fact that F{v) dv < Ce", we deduce: 



]lrf f F\x' -vl'dvdxdt] 
J \Jo Jk" J\v\>s-^ J 



J2 < / F{v)dv 

\J\v\>e-^ 

< Ce" f f f F[\x'\^ + \ipf]dvdxdt 

< Ce" f f( f \x'f + \'pfdxdtdv 

J\v\>e-^ Jo JR" 

< C'e"ll'^lli^((o.oo)x««) / Fiv)dv 

\v\>e ^ 

Incquahty (36) thus yields: 



K{g')[x' ~v] dxdvdt < C£-"/2(e + e") ^ (-(eV^ +£^) 

Jm.'^^ 

and Lemma 4.3 foUows since a G (0, 2). □ 

Proof of Proposition 4-4- We now turn to the proof of Proposition 4.4, which 
we wiU only prove under the stronger assumption 

^(^) = ^fc H^,v) = M^) forallK,|>C, (37) 
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though the result clearly holds in the more general framework (see the proof of 
Proposition 3.2). 

Using (33), we rewrite: 



iy{x, v)F{v) [x'^{x, V, t) — ip(x, t)] dv 



= h+h 
where 



{x, v)F{v) / e- ^0 •^i^+^vs.v) ds^^^ g^^^ ^) r _^ g^^^ ^) _ ^(^^ ^^ 



dz dv 



ii — e 
and 

r _ — ct 

^2 — £ 



\v\<C 



v{x,v)F{v) / e-/o '^(=^+'^''«.")'iV(x+£uz,i;)L(x+euz,t)-(p(x,i) 



\v\>C 



iy{x,v)F{v) / e-i'o ''(^+^'"*'")''V(x+ei;z,z;)L(x+£wz,t)-^(a;,i) 



dz dv. 



We now proceed as in the proof of Proposition 3.2 to show that Ii converges 
to zero, and I2 converges to —C*{if). The only difference with the proof of 
Proposition 3.2 is that the factor 



"(X, V, Z) = e- ^0 ds^(^^ ^ ^^^^ ^) 



p 



is not symmetric with respect to v. 

To prove the convergence of Ii to zero, we thus introduce 

p\x, v, z) = e- ^0 "(^'^^ '^V(x, v) = e-^''(^'^V(x, v), 

and write 

Ii = e^" I v{x,v)F{v) I p'^{x,v,z) Lp{x + evz,t) — ip{x,t) 

J\v\<C Jo 

= e^" / iy{x,v)F{v) / p^{x,v,z) (p(x + evz,t) — ip{x,t) 

J\v\<C 







dz dv 
dz dv 



+e~ 



v{x,v)F{v) I [p'^{x,v, z) — p'^{x,v, z)] if{x + evz,t) — (f{x,t) 

\v\<C 



dz dv 



Thanks to Assumption (Al) and (A2), we have i'{x,—v)F{—v)p'^{x,—v,z) = 
v(x, v)F(y)p'^ (x, V, z), and so proceeding as in (27), we can show that the first 
term is ©(e^""). To bound the second term, wc note that 



and 



\^{x + ewz, t) — Lp{x, i)\ < Ce\v\z 
\p'{x,v,z) -p"{x,v,z)\ < Ce-^'^ielvlz + e\v\z^] 
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We deduce 

< Ce"^-" + Ce-°' [ iy2F{v) [ Ce-''^''[e^\v\^z'^ +e^\v\^z^]dzdv 

J\v\<C Jo 

< Ce^-" 

Where C depends on 1 1£>(/?| |l=° , | |-D^(/3| II-^" , II^k'^Hl'", i^i and 1^2- We deduce 

lim /i = 0, 



It only remains to study the Umit of 12- Using condition (37), we write: 

Jo "'|i'|>C l^'l 

and the change of variable w = ezv yields: 

"'0 J\w\>Cez \W\ L 



dv dz 



lO J\w\>Cez \W, 

Formally, this converges to 



(p{x + w,t) - i^{x,t) 



dw dz 
dw dz. 



PV 



Wc can now proceed as in the proof of Proposition 3.2 to rigorously establish 
this limit. Note that as above, the factor in the integral is not even with respect 
to w. We thus introduce 

—z / i^olx + s{y — x)) ds 
l{x,y,z) = kqVo{x)vo(jj) z"" e Jo 



and split I2 as follows: 

J2 = / / j{x,x + w,z) ip{x + w,t) - ip{x,t) 

Jo J\w\>C6Z 

»1/(C, 

7(x, X + w, z) ip(x + w,t) ~ if{x, t) 



dw 



\N+a 



dz 



J\w\>l 
l/{Ce) 



dw 



\N+a 



dz 



JC£Z<\w\<l 

l/{Ce) 



j{x, X + w, z) ~ 7(2;, x, z) (p{x + u>, t) — ip{x, t) 



dw 



■ dz 



7(x, X, z) ip(x + w,t) — tp(x, t) — Vip{x, t) ■ w 

dw 



1/(Ce) J\w\>Cez 



Cez<\w\<l 

7(a;, X + w,z) ip{x + w,t) — (p{x, t) 



dw 



\N+a 



dz 



\N+a 



dz 
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in order to show that (35) holds uniformly with respect to x and t (all the 
integrals above are defined in the classical sense without need for principal 
value). □ 



A A priori estimates 



We recall here the proof of Lemma 4.1. First, we have the following lemma, 
which summarizes the key properties of the collision operator L: 

Lemma A.l. Assuming that a satisfies Assumption (B2), then the collision 

'2 



operator L is hounded in Lp_i and satisfies: 

I L{f)^dv<-vJ 
where (/) = f{v)dv. 



f-{f)F\^-dv forallfeLl-, (38) 



Proof of Lemma A.l. We adapt the proof of [6, Proposition 1 & 2]. The fact 
that L is bounded in is a simple computation, and the proof is left to the 

reader. To prove the coercitivity inequality (38), we write 

/ L(/)|d« = [ [ a{v,v')fLdvdv'- [ v{v)^-^dv 

JR" ^ JR" JR" ^ JR« -f" 

= / / aiv,v')F'^ldvdv' - [ iy{v)^dv. 

JR" JR" ^ ^ JR« ^ 

Next, we note that the second term in the right hand side can be rewritten 
I v{v)—dv = a{v' ,v)F-—^ dv dv' 

= I I a{v,v')F'Qdvdv', 

JR" JR" ^ 

as well as (using the fact that vF = K{F)) 

( v{v)^ldv = / K{F)^dv 

JR« ^ JR« ^ 

= 11 a{v,v')F' dv dv' . 

JR« JR" ^ 



We deduce 



L{f)^dv = -\ 
^2 



/ / a{v,v')F' 

/R" JR« 


7' 

F' 


/■ 

F 


2 

dv dv' 




/ / a{v,v')F' 
/r« Jr« 


'g' 

F' 


g' 

F 


2 

dv dv' . 


(39) 
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Assumption (B2) then yields 



L{f)jdv < 



2 

2 ./rn 



FF' 



3_ _ g_ 

F' F 



1 2 



dv dv' 



f f F^-2gg' + ^F'dvdv' 



Finally, using the fact that J g{v) dv = 0, we deduce (38). □ 

Using Lemma A.l, we can now prove Lemma 4.1 
Proof of Lemma 4.1. Multiplying (1) by /F, we get: 



2 dt Jr2n ^'^ ^ F 



£} .Ib2N F 



< 



\ff 



.e\2 



R2« 



F 



■ dx dv. 



with = — p^F. We deduce: 
1 f \n 



■ dxdv 



[e) Jo 



F 



< 



■ dx dv ds 



F 



dx dv. 



This inequality shows that is bounded in L°°(0, oo, L'^-i). We also get 



JR2" 



F 



■ dxdvds < C ||/o||l2(^-i) 0{e). 



Finally, Cauchy-Schwarz inequality implies: 



p^fdx 



fdv 



dx 



< I \r?^dv [ Fdvdx= f \f\'ldvdx. 



□ 
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